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Abstract. The most general evolution of the density matrix of a quantum 
system with a finite-dimensional state space is by stochastic maps which take 
a density matrix linearly into the set of density matrices. These dynamical 
stochastic maps form a linear convex set that may be viewed as supermatrices. 
The property of hermiticity of density matrices renders an associated supermatrix 
hermitian and hence diagonalizable. The positivity of the density matrix does not 
make the associated supermatrix positive though. If the map itself is positive, it 
is called completely positive and they have a simple parameterization. This is 
extended to all positive (not completely positive) maps. A general dynamical 
map that does not preserve the norm of the density matrices it acts on can be 
thought of as the contraction of a norm-preserving map of an extended system. 
The reconstruction of such extended dynamics is also given. 
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1. Introduction 

A quantum system with a finitely many dimensional state space may be represented 
by a N x N quantum density matrix p. The density matrix must be of trace class and 
should satisfy the properties of hermiticity and positivity: 

tr(p) = 1 ; jd t = p ; x*p rs x s > 0. (1) 

For a closed system the dynamical evolution of the system is by the action of a unitary 
time-dependent operator. 

p(t 2 ) = U{t u t 2 )p{t 1 )U\t 1 ,t 2 ) (2) 

where 

ta 



U{tx,t 2 ) = T jcxp ^-t jf 2 H{t')dt^j | 



The evolution is linear. But if we have an open system and if we are considering 
the back reaction of the environment on the system then the dynamics cannot be by 
unitary evolution but by a more general linear evolutionpP: 

p(t 2 ) = A{t x ,t 2 )p{tx)- (3) 

The linearity of A(t\,t 2 ) follows from the linearity of quantum mechanics and for 
this reason we do not consider more complicated forms of maps on density matrices. 
The superoperator A can be written as a supermatrix and the transformation can be 
written as 

Prs * A rs - r i s ' p r i s ' (Ap^ rs . 

In the equation given above, the elements of the density matrix has been suitably 
regrouped into a column vector so that A can be in the form of a N 2 x N 2 matrix. 
The constraints Q on the density matrix impose restrictions on A. It is instructive 
to first recast A into another dynamical matrix 13 |S] B such that 

Ar s;r' s' (t*) B rr '- S ' s (t^. 

In the form B the supermatrix has to satisfy the following relations: 

Ks^r'it) = B rr , ;s , s (t) (Hermiticity) (4) 

B n r'-s'n = b"r's' (Normalization) (5) 

x* r y r >B rrl . s , s x s y* s , > (Positivity). (6) 

In terms of the pairs of indices rr' and s's, B is a hermitian matrix which gives 
non-negative expectation values for factorizable supervectors. 

v) Bu > if u rs = x r y*. 

It is not necessary that B > for maintaining the positivity of the density matrices 
under dynamical evolution even though it is a sufficient condition. If B > we will 
call the map "completely positive" |2j • The terminology is slightly confusing unless one 
keeps in mind that positivity of the map is a statement about its action on density 
matrices while complete positivity can be regarded as a statement about the map itself, 
in addition to saying something about its action. For instance the action of taking the 
transpose of a density matrix: p — > p T is a positive map but not a completely positive 
map. Complete positivity ensures that the evolution of the system of interest, S due 
to the action of the map is extensible in a trivial way to a physical evolution in a larger 
system 5" <g> S. Such an extension is especially useful when S is entangled to S' . 
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Completely positive maps are interesting in the context of quantum information 
theory and quantum computing in that it provides a description of the types of 
dynamical evolution of a composite quantum system that can change the degree of 
entanglement between parts of the system. In other words we can deal with the open 
evolution of quantum systems rather than unitary isolated evolution thereby getting 
a handle on the difficult problem of the influence of the environment on a potential 
quantum computer. Extensive literature exists on the subject; we refer the reader to 

IZ] and the references therein. Completely positive maps have also been used in 
certain cosmological models and in attempts to construct a quantum theory of gravity 
in the context of Quantum Causal Histories joj. 

In this paper we look at maps which are positive but not completely positive. 
Not completely positive maps do not represent dynamical evolution in the sense that 
completely positive maps do. This point is discussed in sections 01 and 0] 

The density matrices on which the dynamical maps act form a convex compact 
set. Since positive maps transform such a set into itself, they themselves form a convex 
set (the set of completely positive maps also form a convex set lying inside the set 
of positive maps). The convexity property means that any linear combination of of 
positive maps with nonnegative coefficients which sum to unity is also another valid 
map. For instance the positive map B defined by 

B = Yl k ( n ) B n ; k(ri) > ; fc W = 1 

n n 

is also a positive map if B n are positive maps. It can also be shown that 

1 < n < N 2 where N is the dimensionality of the density matrices that the map acts 
on. Out of the convex set of positive maps we can pick out those maps which cannot be 
written as a sum of other maps. Such maps are called extremal. In the discussion that 
follows we talk only about extremal maps and is assumed to denote such maps. 
Any generic map can be constructed out of the extremal maps as a linear sum through 
the specification of at most N 2 extra parameters. We look for parameterizations of 
extremal maps. This, in conjunction with the N 2 extra parameters which add up to 
unity, will then suffice to parameterize a generic dynamical map acting on density 
matrices. Note that such a characterization of the extremal maps is usually carried 
out with the assumption that the additional restriction in equation (5) that the maps 
be trace preserving also holds. This restricts the maps being considered to a bounded 
convex set rather than to a convex cone. 

Since B is hermitian it follows that it has an eigenvector decomposition 

* (a ?C&>* = CMQ\ (7) 



where M is a diagonal matrix with eigenvalues fi a and Q?/ are the normalized 
eigenvectors. For a completely positive map all the \x a are nonnegative but this is 
not true for all positive maps. If all the /x Q are nonnegative, we can absorb them into 
the eigenvectors by defining 

r {a) - nr r (Q) 

rr' V r la Srr' 

So for a completely positive map 2 
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with the trace condition 

a 

Parameterization of extremal completely positive maps can be found in reference |17| . 
Note that if the completely positive map is not extremal, its action can be written as 

p^Y. k ^ C n ] P C( n )] - ( 9 ) 



2. Not Completely Positive Maps 

The action of an extremal map which is not completely positive on a density matrix 
can be expressed as 

m n 

p^Y. c[a) p c(a)] -Y, Di0) p Dm ( 10 ) 

a=l 0=1 

with 

where v are the negative eigenvalues of the map and vffj are the corresponding 
eigenvectors. 

It will turn out that the number of positive eigenvalues m have to be greater than 
or equal to the number of negative eigenvalues n. The trace condition now becomes 

m n 

j2 crwtcw - ]T dWdW = i NxN . (ii) 

a=l 0=1 

The positivity condition yields the general result 

m n 

J2 C^u a ulC^ - J2 D^upu\D^ > (12) 

a=l 0=1 

Taking the trace on both sides of equation Ijlll) we obtain: 

m n 

J2 tr [cWtc-W] - tr [d^dW] = N (13) 



0=1 



Define 



C ^C {a) = J > (14) 

a=l 
n 

Y^D^D^ =K>Q. (15) 
0=1 

Since C (a ^C^ is Hermitian so is J. We can therefore perform a unitary 
transformation {/ that diagonalizes J. By the trace condition the same unitary 
transformation automatically diagonalizes K . Then equation 1)11(1 becomes 

j- k = i 

where J = [/ JC/ t and A" = [/ifl^ 
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Let the eigenvalues of J be jf and those of K be fc 2 (0 < i < N). Since 
J = diag(j 2 , j\ . . . jfj) and K = diag(fc 2 , k\. . . k 2 N ) we have the relation 

j} = kf + l 

from the trace condition. Since kf > the number of eigenvalues ji must be greater 
than or equal to the number of ki, i.e. to > n as mentioned earlier. 
Define tpi so that 



Now define 



ji = cosh (fi , 

C (Q) = [coshtp] M (a) 
# = [sinh^] NW 



ki = sinh ipi . 



(16) 
(17) 



where we have extracted the matrices [coshyi] and [sinh<ys] from the matrices and 
£)(0) respectively. It follows from equations Ijl4|l and (|15|l that 



A/ (Q)t M (a) = 1 

a=l 



(18) 



J2 nVNnW = 1. 
0=1 



(19) 



Parameterizing the matrices M and N are already known. In the cases where 
sinh<^ = 0, for all tp, we have only a smaller set of matrices to parameterize and 
this is identical to the case of having a completely positive map. Here we assume 
that sinh tp and see how many parameters we need to write the map in the most 
general case (up to a unitary transformation). 

Since the matrices M^M are hermitian we first choose a unitary transformation 
W\ that can diagonalize M^M^: 



W\M^M {1) W X = 



o2 fl(l) 



\ 



cos w 2 



V o 



cos 2 ^ J 



We can make further simplifications on MWtMW by noting the following P^: 
We have at our disposal now m matrices such that 



E 

a=l 



(20) 



out of which the first term, WIM^ 1 ^ W\ has been diagonalized. We can now 
make the following transformations on the first two terms of the sum in l|2U|l without 
changing the sum itself or the eigenvalues of the map: 

WlM^M^Wi -> 



( cos 2 df - 1 



W}M^M^Wi 







\ 



o / 
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WlM^M^Wx 



WlM^M^Wi 



COS" 



V 







o \ 



o / 



Such transformations which may be done on the set of matrices {M} (or on {iV})do 
not change the sum in l|18f) (or in (|19l) ). The map which is made up of the sums 
in (|18() and (|19f) hence remains unchanged and so these transformations allow us to 
fix cos 9^ = 1. This freedom corresponds to performing orthogonal transformations 
(rotations) in the space containing the matrices M^ a \ This in turn can be interpreted 
as an orthogonal transformation on the density matrix p on which the map acts. An 
alternate way of looking at the transformation is to note that the sum of two matrices, 
A + B can always be written as the sum of two other matrices say, for instance, C + D 
where C = (A + W)/2 and D = (A — B)/2. Once we have set cos 9^ = 1, we 
parameterize using N — 1 angles ; 2 < i < N. 

Applying the transformation W% to both sides of equation l|18|) . exploiting the 
extra freedom mentioned above and after doing some re-labeling, we obtain 

/ o \ 



M[ a) ^M[ a) 



sin 2 0« 







9 (1) 



(21) 



Where M[ a) = W^M^Vl for 2 < a < m. V\ is another unitary matrix that reduces 
M" to the form that we want. We can now focus on the set of (N — 1) x (N — 1) 
matrices defined as 



M 



(a) 



V 



sin 



i(i) 



\ 



sin 



3(1) 







sm 



i(l) 

7 N 



M[ a) 



2 < a < N 



where we have dropped the first row and column of on the right hand side of 

the equation and also extracted the factor containing sin 6^ from it. We assume that 

(i) 

none of the sin 9\ are zero since we are interested in computing the maximum number 
of parameters required for describing a generic extremal map that is not completely 
positive. 

From equation (|21|) it follows that 



a=2 



(2) 

the matrix M{ in the first term of this sum can be parameterized using N — 
2 parameters using exactly the same procedure as before. Using a unitary 
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transformation W% and a further orthogonal transformation (if needed) we can 

(2) 

transform M{ to the following form as before 

/ 1 



Afj 



(2) 



cos 



,(2) 







\ 



Repeating this procedure m times, we parameterize all the matrices C^ a \ can 
also be parameterized in the same fashion. The total number of parameters needed 
can be computed as follows. There are N angles tpi. To parameterize the matrices 
M (q) we need (N - 1) + (N - 2) + . . . (N - m) parameters and for N {m we need 
(N — 1) + (N — 2) + . . . (N — n) parameters. So in total we need 

m(m — 1) + n{n — 1) 



N 2 - 



(22) 



parameters. 

Note that the matrices CM and are determined only up to m + n unitary 
matrices according to 

(j( a ) > (j( a )jj( a ) 

which leave equation l|13l) unchanged. We can see this also in the manner we defined 
M[ a) = W\M^ V\ where we had to introduce the arbitrary unitary matrix V\ . 



3. Dynamical Maps as Contractions 

A straightforward way of generating positive maps is to consider the unitary evolution 
of two systems coupled to each other. Let S be the system of of interest and R the 
second system. Let the dimensionality of S be d and that of R be N. R can be treated 
as a 'reservoir' with which S is interacting. Dynamical maps representing the time 
evolution of S can then be thought of as contractions on the unitary evolution of the 
combined system. If we choose a direct product density matrix as the initial state 
then the dynamics of the coupled system is given by 

11 = psxtr — > Vp s X T R V^ 

where V is a unitary matrix in the direct product space Tls x Using the index 
notation employed in the previous discussion 

Prs x Tab > V ra - r ' a ' p r ' s ' T a i\,i V sb . s i b , . (23) 

The evolution of the system S is extracted using the partial trace operation which is 
a contraction. 

Prs > tr^ ( Vra;r'a' Pr' s' T~a'b' ^sb-.s'b' ) — Vrn\r' a' Pr' s' Ta'b' Vsn-s' b' ' 

For simplicity we assume that r can be made diagonal by a suitable unitary 
transformation in Tin with eigenvalues t(1), t(2) . . . r{n). Then the map on S is 

Prs — ► ^2 V rr >(n,v)p r > s >T(v)V*< s (n,v). (24) 

v.n,r' ,s' 
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Here the operator V has been rewritten in a manner suggestive of the form of a 
completely positive map that is not extremal, given in @ . i.e. 

v a 

To get the form of the map in 10 it is sufficient that the dimensionality N of the 
reservoir to be the same as that of the system, i.e. N — d. With this restriction t(v) 
could correspond to a mixed state. If we further restrict r to correspond to a pure 
state so that it has only one eigenvalue then the map is extremal and reduces to the 
standard form (jHJ 

p rs = £ vMp^vW* * c {a) P c^ 

n a 

in which a runs over 1 < a < d. In other words extremal completely positive maps 
are contractions of unitary evolution in a space in which the system is coupled to a 
reservoir, of the same number of dimensions of the system, whose initial state is a 
pure projection. 

We note here that all these maps are completely positive maps (not necessarily 
extremal), if the dimension of the reservoir is made bigger than or equal to d 2 then 
any map on S (not necessarily extremal) can be expressed as a contraction of the 
unitary evolution in a space in which the system is coupled to a reservoir whose initial 
state is a pure projection. 

We can also carry out the inverse construction where we start with a completely 
positive map and view it as a unitary transformation on a larger system. Given an 
extremal completely positive map of the form © , we can construct a unitary matrix 
V in mn dimensions with: 

\r — r<( Q ) 

Vra-,r'l — ^ rr' ' 

The conditions on C 1 -"' are transcribed into 

which is necessary for V to be a unitary matrix. The ambiguity in constructing the 
other elements of V, where the last index is not equal to 1, does not affect the map. 
In short, V can be constructed in such a fashion that it corresponds to any given 
dynamical map on the system along with a particular choice of the states r of the 
reservoir. 

4. Not Completely Positive Maps as Contractions 

What about not completely positive maps? To obtain such a map as a contraction 
we have to generalize the auxiliary space TIr to be a space with an indefinite metric 
and V to be a pseudo unitary operator in the mn dimensional space. Positivity of the 
map is guaranteed if the generalized density matrix of the extended system is initially 
entirely within the convex set of positive metric states of the mn dimensional space. 
The sum over the index n in l|24[) goes over both positive and negative metric terms; 
but the resultant density matrix is nonnegativc. In general the not completely positive 
maps have to be viewed as contractions of the evolution of unphysical systems. For 
example the complex conjugation of the density matrix p of the system is a positive 
map which is not completely positive. Physically it has the meaning of time reversal 
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of the system. However if we view it as a contraction of the evolution of two coupled 
systems, then this corresponds to time reversal of only of the systems which is rather 
meaningless. 

We can invert this derivation to realize the most general not completely positive 
map as the contraction of a larger evolution in an indefinite metric space for the 
reservoir. To make the map extremal we further restrict the density matrix of the 
reservoir r to have a single eigenvector (with positive metric) with eigenvalue unity. 
Since such reservoirs are somewhat artificial, we have to consider this reconstruction 
as a purely formal device. 

5. Summary 

We have studied linear dynamical maps which take the set of density matrices into the 
set of density matrices. These maps form a convex set which is also compact in the 
case of completely positive maps. It has already been shown in j!7j that a completely 
positive extremal map contain at most N terms, requiring a total of N(N — l)/2 
terms to parameterize each of the terms up to a set of N x N unitary matrices. The 
completely positive maps can be viewed as the contraction of unitary evolution in an 
extended space. Extremal maps correspond to the case where the auxiliary system in 
the extended space is a pure projection. Conversely we can reconstruct the unitary 
evolution of the expanded system from the map itself. 

These considerations are extended in this paper to positive but not completely 
positive maps. The extremal maps again have at most N terms. The number 
of parameters required to describe each one of these terms up to a set of unitary 
N x N transformations is given in equation We can obtain these maps also 

as contraction of an extended system. But here the extended system has a pseudo 
unitary evolution matrix. It is also possible to obtain this pseudo unitary evolution 
starting from the maps. 

These results generalize the results obtained two decades ago by Gorini and 
Sudarshan^J for 2 x 2 matrices. 

Needless to say, however complicated the dynamical processes leading to the linear 
stochastic evolution that is represented by the dynamical map, we see that the same 
dynamics obtains when we couple the system to a reservoir having dimension N 2 . In 
the case of an extremal map it suffices to couple the system to a reservoir of dimension 
N. 

In this paper we have dealt only with dynamical maps and not with the continuous 
semigroup of evolution. This study was carried out by Kossakowski 16 and followed by 
others |13U14| . In these, while the semigroup generators are parameterized no attempt 
is made to embed them in a larger system. Since the Zeno effect^! operates for very 
small time intervals, care must be taken in generating a semigroup from the dynamics 
of an extended system. We hope to examine this question in the near future. 
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